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O ■ The environments of supernova explosion are often inhomogeneous and there 

may be jumps in their density structure. We have developed a semi-analytic 
model under the thin-shell approximation for supernova remnants that evolve 
crossing a density jump in the ambient medium. The generic evolutionary rela- 
tions are presented for the blast wave after impacting on a cavity wall, which may 
be produced by the energetic stellar wind from the supernova progenitor. The 
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relations can also be extended to the case that the blast waves break out from a 
dense cloud if different density contrast is used. This model is applied to N132D, 
a well-known cavity-born supernova remnant whose evolution has not yet been 
quantitively estimated in a cavity scenario due to lack of model formulae, and 
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self-consistent physical parameters are obtained. 

Subject headings: ISM: bubbles — ISM: kinematics and dynamics — shock 
waves — supernova remnants — supernova remnants: individual: N132D 



^ . 1. Introduction 



The environments of supernova explosion are often inhomogeneous and there may be 
jumps in their density structure. A massive star can excavate a low density cavity with its 
energetic stellar wind and ionizing radiation in the circumstellar space before it explodes. 
It is also possible that supernova may explode in a dense cloud and subsequently break 
out from the cloud into a low density region. These supernova remnants (SNRs) cannot be 
simply treated with the canonical evolutionary laws such as the Sedov relation (Sedov 1959). 
In the former case, a blast wave collides with a cavity wall after it expands effortlessly into 
the cavity. A similar situation has been suggested to have occurred in various supernova 
remnants such as N132D (Hughes 1987). In the latter case, the blast wave propagates at a 
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higher velocity after it breaks out of the cloud, as is suggested for 3C 391 (Reynolds & Moffett 
1993). It would be of interest to develop an analytic model showing how the evolution of 
such kind of supernova remnants deviates from the canonical solution. 

In a study of SNR evolution in a dense molecular cloud, Wheeler et al. (1980) briefly 
discussed the effect of a pre-SN cavity, suggesting that the remnant will evolve rapidly 
into the radiative phase after it encounters the cavity wall. Chevalier and Liang (1989) 
presented a self-similar model for both the structure of the shocked powerlaw ejecta and 
the shock propagation into the circumstellar shell. However, the ejecta structure becomes 
unimportant when the mass swept up by the blast wave significantly exceeds the mass of the 
ejecta. Investigation of supernova blast waves crossing a density jump has been carried out 
by hydrodynamic simulation. Tenorio-Tagle et al. (1990, 1991) made numerical simulations 
of the supernova shock interacting with a wind-driven shell (see also a review by Franco 
et al. 1991 and references therein). Tenorio-Tagle, Bodenheimer, & Yorke (1985) have also 
presented a numerical simulation of the evolution of a remnant resulting from supernova 
explosions in or near molecular clouds. 

An analytical model would not only complement the numerical simulation, but would 
allow for exploring a large parameter space and more directly provide physical insights into 
the evolution of this sort of SNRs. It has been shown by many authors (e.g., McCray 1987, 
Blinnikov, Imshennik, & Utrobin 1982, and others) that a thin-shell model (Kompaneets 
1960), though approximate, is convenient and valid in revealing the general evolution of 
the supernova blast wave. The approach was applied in a numerical algorithm to calculate 
the evolution of SNRs in a density gradient (see Wang et al. 1992). This method has also 
been applied analytically to other explosive events, e.g., a cosmic ray blast wave (Morfill 
& Drury 1981). Therefore we use the same approach to obtain an analytic/semi- analytic 
approximation to seek the basic evolutionary rule of the blast wave after an impact on a 
cavity wall or a breakout from a cloud. 

2. Thin-shell Model of Impact on a Wall 

For the ambient environment of supernova, we assume a density structure of medium 
p = Pi = lAnim H for r < r j and p = p Q = lAn Q m H for r > r j (where rii and n are the number 
densities of H atoms inside and outside the density jump, respectively). In this scenario, the 
ejecta or blast wave is assumed to expand easily until it hits upon the cavity wall at r- y By 
the time the ejecta or blast wave arrives at rj, the SNR could be in the Sedov phase or even 
still in the free expansion phase, depending on the parameters rj and pi. Denote the density 
ratio by j3 = (pi + p C j)/p , where p e j is the average density of the ejecta mass over the cavity 
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within rj, so here we discuss the case (3 < 1. The swept-up material is assumed to be in a 
thin shell and the mass in a conical section of solid angle Ail is given by 

AM B = ±Mlpo[r*-rf(l-j3)]. (1) 

On the assumption that most of the explosion energy is thermalized by the time when the 
blastwave impacts the wall or soon after that, the thermal pressure in a spherical remnant 
is given by 

= ( 7 -l)£ = _E_ , . 

P Anr 3 s /3 27rr s 3 ' 1 ' 

where E is the total energy and the adiabatic index 7 = 5/3. The momentum equation for 
the conical section is 

(3) 



-(AM s v s )=pr 2 AQ, 



From the equations (1), (2), and (3), we have 
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2.1. Adiabatic phase 

Before the radiative phase, the supernova remnant is assumed to be adiabatic and the 
energy to remain unchanged: E rs E . Then eq.(4) has a solution 

«l=( ^rV= ( —) ^"'"'^'f. (5) 



dtj \npoJ [ r 3_ r 3(i_/3)]- 

where integral constant C must be determined according to the initial condition at r s = r- y 
The problem can be divided into two cases presented below. 

In the first case, the supernova remnant has already been in the Sedov phase by the 
time tj when r s — rj. The Sedov evolution law under the thin shell approximation is given 
by r s = (£E / pi) 1/5 t 2/5 where £ = 25 /4n (e.g. McCray 1987), so the initial velocity at r- s is 
v 2 = E /(npirf). The blast wave velocity is thus 

^ ^T ^s 3 -r J 3 (l-^)[l-31n(r J /r s )] 

dt \jn Po fj-fjti-p) (6) 

namely 
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V2 ^ (A) (Aj), (7) 
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where Aj = Tj/r s and 



/ 1 _ A?(1 _ /J)(1 _31 IlA ) 



In the second case, the remnant is still in the free expansion stage when arriving at rj. 
In this case, pi < p c y, and, for simplicity, /3 ~ is assumed. For /3 — 0, the above initial 
condition cannot be used, and instead, the convergence of eq.(5) at r s = r- } is entailed, so 
C = — rf + 3r?lnrj. Thus, interestingly, the velocity is still described by eqs.(6)-(8) (but 
with (3 = 0). 



Eq.(6) gives rise to 



f rs r 3 -rf(l-p) 



or 



where 



t~t 3 = J^ ^ = dr (9) 

V A> J rj Jr 3 - rf(l - 0)[1 - 3 ln(rj/r)] 

[^ (A) (A J )] 2/5 r s = (t-t,) 2 / 5 , (10) 



S/2 



^, (A) (Aj)= / . A ^ d\= / dA (11) 

A ,/A 3 -A3(l- / 3)[l-31n(A j /A)] J\ ^ (Aj/A) 



j 

with A = r/r s . The above integration (at rj 7^ 0) can only be performed numerically, as 
plotted in Fig.l. 

Note that, for the special case rj = Aj = tj — 0, F^ A) = 1 and F r (A) = 2/5, so we 
immediately come from eq.(10) back to the Sedov solution r s = (£E / ' p ) l ^t 2 ^ where £ = 
25/47T under the thin-shell approximation (e.g. McCray 1987). 

As shown in Fig.l, both functions i^ A) (Aj) and F r (A) (Aj) vary with parameter (3. For 
(3 7^ 0, function i^ A) (Aj) starts from /3 -1 / 2 , as a limit value found for eq.(8) at Aj = 1 or 
r s = Tj. The remnant is decelerated rapidly after it enters the cavity wall, as indicated by 
the ensuing rapid decrease of i^ A) (Aj). Just when r s surpasses r- v the smaller (3 is, the more 
drastically F^ A) drops and the more quickly F^ A) approaches the asymptotic curve for (3 — 0. 
In addition, these asymptotic values are of order unity, implying a crude factor of velocity 
decrease ~ (3~ l l 2 or a ram pressure (pv 2 ) roughly recovered to that before the impact (see 
eq.[7]). 

In the limit case (3 — 0, immediately after the remnant strikes the cavity wall [1 — Aj = 
(r s — rj)/r s <C 1], the behavior of F^ A) can be described by the Taylor series 



l + ^l-Aj) 
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and the shock velocity at r s = rj (i.e. Aj = 1) is 

Va (r i )= f-g^V /2 = 1,4 x 10 3 f s f° V /2 f^-r V Y^ km s- 1 , (13) 
v j; V 2 wf/ \10 51 ergs/ Vlcm- 3 / V 5 P C / 

a factor lower than the velocity v s = (E /itp r?) 1 / 2 the shock would otherwise have in 
the Sedov solution. Still in the case (3 — 0, the linear relation [eq.(12)] is found to be a good 
approximation of F^ A) for 1 > Aj ^0.2 (see Fig.l). 

At the other extremity Aj <C 1, F,J A) (Aj) approaches to 1 according to 

^ A) (A j )^l + i(l-/?)(l + 31nA j )A J 3 , (14) 

which allows for any values of (3. 

The numerical solution of function F r (A) (Aj) in eq.(ll) shows that it starts from at 
Aj = 1 and approaches the other extremity Aj = in similar trends for various values of (3 
(Fig.l). For 1 — Aj <C 1, in the case (3 ^ 0, its behavior can be approximated by a series: 

F^{\) » ^= [(-9 + 13/?) + 50(1 - (3)\ - 105(1 - (3)\] + 4(16 - ll(3)\f] , (15) 

which is obtained by expanding function i^ A) (Aj/A) in eq.(ll) to the second order of the 
small quantity (1 — Aj/ A); in the case (3 = 0, however, 

)(Aj) « (3 + 2A, + 3Af - 8Af) . (16) 

This series expansion is found valid in a broad range 1 > Aj k, 0.3 (Fig.l). 

At the extremity Aj <C 1, function F r (A) (Aj) converges to 0.4 conforming to 

F r ^(Aj) * \ (l - Af) + 5(1 - (3) (Af - A 3 ) + 3(1 - /?)Af In Aj, (17) 

which is valid for any value of /3. 

For (3 ^ 0.1, numerical integration (Fig.l) shows that F r (A) (Aj) as a whole does not vary 
significantly with various (3. If we take (3 = 1 for simplicity, it corresponds to the uniform 
case; eq.(8) yields F^ A) = 1 and eq.(ll) yields 

^'(AjH^l-Af) (18) 

(see the long dashed line in Fig.l). If this expression is adopted as a crude approximation 
of F r (A) (Aj) for (3 ^ 0.1, eq.(10) then gives a Sedov-like evolution 

'-^Jiki'T-ry (») 
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where £ = 25/47T but can be replaced with the canonical value 2.026. 

The adiabatic phase will come to an end when the gas temperature drops to T c = 
6 x 10 5 K (Blinnikov et al. 1982). With the relation kT s = (3/16)/im H v 2 (where the mean 
atomic weight \i = 0.61), Eq.(7) with p = lAn m H gives the location r c of the phase 
transition: 

' -(^) 3 (J^:)(^^r (20) 



[F^ir-Jr,)] 2 Vl0 51 erg S ; Vl0 2 cm- 3 

If (3 <C 1, the linear approximation (12) is adopted and the above formula is simplified as 

- 2 r-j « (3.8 pc) 3 E 5 i (n /10 2 cm" 3 ) _1 (21) 



r r .r\ 



where E§\ = E / (10 51 ergs). For (3 <C 1, the sudden deceleration of the remnant on impacting 
the cavity wall makes it possible that the shocked gas temperature T s (rj) has already gotten 
below T c and the radiative phase has begun before the velocity falls on the asymptotic curve 
(for f3 — > 0, see Fig.l). In this case, the adiabatic phase lasts so little time after the impact 
that it essentially may be ignored. Therefore it is possible that a previously freely expanding 
remnant directly enters the radiative phase after it strikes the wall and the adiabatic phase 
is abortive if rj > 3.8 pc El( 3 (n /10 2 cm" 3 )- 1 / 3 . 



2.2. Radiative phase 



When the shell is cooled down by radiation, it is assumed to be driven by the thermal 
pressure of the still hot gas interior to the shell. The internal hot gas suffers from the energy 
loss by the expansion: 



dE_ 
~dt 



—4:7ir 2 v s p. 



Eq.(22), inserted with eq.(2), yields 



E = E {r c /rf 



Now, eq.(4) becomes 



o d 
r — 

s dt 



SE r 2 c 
2np u 



(22) 



(23) 



(24) 



with a solution 

dr R l3E^ c \/^ - r c) + ^f(^ 2 - r- 2 ){l - [3)/2 + [r c - (rf/r?)(l - (3) (1 - 3 ln^/rj)] /3 



dt 



rl-rf(l-P) 



(25) 
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In this solution, we have used eq. (6) as the initial condition of the shock velocity at the phase 
transition point r s = r c and t = t c . 



Setting A c = r c /r s , r\ = r c /rj, and fj, — (1 — j3)/rj 3 , eq.(25) is rewritten as 

„2\ 1/2 



where 



F(R)(Ac) = V(l - A c ) + (/i/2)A c (Ag - 1) + (A c /3)[1 - Ml + 31nq)] 



1 - M c 3 



(26) 



(27) 



Eq.(25) gives 



3£ r 2 



x 



x 



~ 3 -rf(l-/3) 



or 



_ rc ) + r 3( r -2 _ r -2) (1 _ py 2 + f c _ ( r 3/ r 2)(i _ J) (i _ 3 ln(rj/r c ))] /3 



[Fr(Ac)f /7 r s =f^y /7 (t-tc) 2/7 , 



dr 
(28) 

(29) 



where 



F r (R) (A c ) = 



A 3 - /iA; 



a c V(A - A c ) + (/i/2)A c (A 2 /A 2 - 1) + (A c /3)[1 - //(l + 3 In??)] 

1 A 5/2 



rfA 



a c ^ r) (a c /a) 



dX. 



(30) 



The numerical values of functions i^ R) (A c ) and F r (R) (A c ) are plotted in Figs.(2)-(4), 
where it is seen that both functions depend on the parameters (3 and r\. For most cases 
of (rj,/3) combination, Ff» (A c ) starts from [1 - /i(l + 3 lnr])] 1 / 2 /[ v / 3(l - //)] at A c = 1 (i.e. 
r 's = r c ); but for (r),/3) = (1,0), it starts from and the Taylor series 



^ (R) (A C 



l + ^(l-A c ) + |(l-A c ) 2 



(for 1-A C «1) 



(31) 



can acts as a good approximation in the range 1 > A c k, 0.5 (Fig.2). In the latter case, when 
the radiative phase begins at r s = r c = rj (Aj = A c = 1), eq.(26) turns back to eq.(13). For 
any (77, /3) combination, F^ (R) (A C ) approaches 1 at the other end, A c = 0, according to 



Fr(A c )~l--[4 + //(5 + 61n ?7 )]A c . 



(32) 
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It is noted that, at the limit A c — > (or r s > r c ), one has F^ R) = 1 and F r (R) = 2/7, 
and then eq.(29) is identical to the canonical solution for SNRs in radiative phase: r s = 
( U7E r 2 c /A7ip ) 1 l 1 t 2 l 1 (e.g. McCray 1987; Blinnikov et al. 1982). 

While we could not find an analytic integration for eq.(30), we seek for a series approx- 
imation for A c <SC 1: 

^- (R) ( A c) ~ I 1440 + 168A C (4 + 5/i) + 35A 2 (4 + 5p) 2 - A 7 / 2 (2672 + 2240/i + 875/i 2 ) 
5040 

-84A c/ u [-12 - 5A C (4 + 5/i) + A^ /2 (32 + 25/j)] In rj -1260A 2 (-1 + X 3 J 2 )/i 2 In 2 ??} (33) 

Eq.(33) is obtained by expanding the integrand in eq.(30) to the second order of small 
quantity A c /A and, for the case f] ^ 1.1, is found to be very close to the numerical values in 
the entire range (see Figs. 3 and 4). In the case r\ = 1 and ,8^0, similarly we get a series 
approximation for 1 — A c <C 1: 

F™(\ c ) « [(3 (135 - 434A C + 455A 2 - 156A 7 / 2 ) + 7 (-15 + 54A C - 55A 2 + 16A C 7/2 )] , 

(34) 

and in the case (77, f3) = (1, 0) (1 - A c < 1), 

(Ac) « H (27 + 42A C + 35A 2 ) - |y| A c 7 / 2 . (35) 
Eq.(35) is a very good approximation for (77, (3) = (1, 0) in the range 1 > A c ^ 0.3 (Fig. 2). 



3. Thin-shell Model of Breakout of a Cloud 

We assume a configuration of environment density similar to that in the former case, 
except that pi > p or f3 > 1. The material in the cloud is dense enough that we can assume 
the supernova remnant has been in the adiabatic (Sedov) phase by the time of breakout. 
Moreover, because the density p outside the cloud is so low, we need not consider the 
radiative phase in this scenario, as it is unlikely the SNR would yet have entered this phase. 
Therefore most of the evolutionary solutions for the adiabatic phase established in the former 
case are compatible with this case, given a replacement of f3 > 1. Functions i^ A) (Aj) and 
F r (A) (Aj) for 13 > 1 are plotted in Fig. 5. 
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4. Applications 
4.1. Applicability 

The model developed above can, strictly speaking, only be applied to an SNR of spherical 
symmetry. In reality, an SNR is usually not symmetric, which could be caused by the 
motion of the progenitor for instance. While a good fraction of stars show a considerable 
proper motion, there are still objects with a low proper motion. In practice, as long as the 
asymmetry is reasonably small, our model should still provide a useful approximation. As 
a general condition, the dynamical response (characterized by a sound traveling timescale) 
due to the non-spherical symmetry of the density jump surface should not be rapid enough 
compared with the blastwave propagation after crossing the jump surface. 

Furthermore, we have assumed that the density jump happens only once. However, this 
may not be true for an SNR evolving in a stellar wind bubble, for example. In this case, a 
dense shell, surrounding a low density bubble, may be approximated as two density jumps. 
In principle, such multiple density jumps can be incorporated into an analytic approach. But 
this is a bit too complicated to be included in the present work. Here, we concentrate on 
exploring the effect of one density jump on the evolution of an SNR, which is still evolving 
in the shell. This particular situation is important, because of the expected strong density 
jump effect on the X-ray emission of the SNRs. 

In the range of applicability discussed here, the solutions derived in section 2 can con- 
veniently used in estimating the parameters of relevant supernova remnants. For example, 
X-ray observations provide the X-ray luminosity or the volume emission measure of the X- 
ray emitting gas as well as the gas temperature that is related to the velocity of the blast 
wave. Note that the temperature measured in a plasma model fit to the overall spectrum 
of an SNR only represents a characteristic mean temperature of the shocked gas, even for 
the Sedov blast wave. In principle, one can derive a relation between this mean temperature 
and the shock velocity. This problem may also be resolved to some degree with the spatially 
resolved spectroscopic data, which allow for the extraction of the spectrum of the freshly- 
shocked gas. By use of eq.(7) (or eq.[26]) and eq.(10) (or eq.[29]) together with a relation 
between the volume emission measure and the density of the ambient medium, any three 
among the five parameters, the density n , explosion energy E, the age of the remnant t, 
the density contrast f3, and the fractional radius of the density jump Aj, can be derived from 
the other two, given the three known parameters, the radius remnant r s , the shell velocity 
v s , and the volume emission measure EM. According to this recipes, we apply the model 
developed above to SNR N132D that are conjectured to have expanded from a stellar wind 
cavity into a region of dense medium and derive relevant physical parameters. 
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4.2. SNR N132D 

N132D, an SNR in the Large Magellanic Cloud, lies near the northern boundary of an 
associated molecular cloud (Banas et al. 1997). The inconsistency between the kinematic age 
and the Sedov dynamical age led to the suggestion that the supernova explosion occurred 
in a low-density cavity evacuated by the stellar wind and ionizing radiation of a high-mass 
progenitor (Hughes 1987; Morse et al. 1996; Hughes, Hayashi, & Koyama 1998). Morse et 
al. (1995) derived a kinematic age of 3150 ± 200 yr from the measurements of the velocities 
of fast moving oxygen-rich filaments. Morse et al. (1996) studied the optical photoionization 
precursor based on the [OIII]A5007 surface brightness and derive a preshock density of n ~ 
3 cm -3 and a blast shock velocity v s ~ 800 km s -1 . They find that these results would be 
trapped in unreconcilable contradictions (either the inferred explosion energy is unusually 
high ( ?Z 10 52 ergs), or the inferred age is much larger than the kinematic age), unless the blast 
wave has been traveling much faster in the past, probably by expanding into a low density 
cavity. However, quantitive estimates of the evolution of the SNR expanding from the cavity 
into the cavity wall have not yet been carried out because of lack of proper evolutional 
formulae. 

N132D is an oxygen-rich SNR but no evidence was found by the Hubble Space Telescope 
(HST) for the presence oxygen-burning products. This fact was taken by Blair et al. (2000) 
to suggest a Wolf-Rayet (WR) star which ended at a Type lb SN explosion. However, X- 
ray data from XMM-Newton (Behar et al. 2001; Aschenbach 2002) do show the presence 
of significant proportions of Si and Fe. If these elements are actually present in the ejecta, 
rather than the surrounding ISM, this might tend to contradict the scenario of a Type lb 
explosion. Earlier X-ray findings by Hwang et al. (1993) using the Focal Plane Crystal 
Spectrometer on Einstein also found strong Fe lines, which the authors used to support the 
scenario of a massive progenitor (with masses over 2OM favored). Therefore we assume 
that the core collapse SN explosion occurred in a stellar wind cavity blown by a massive 
progenitor, possibly a WR star. 

The recent XMM-Newton X-ray observation of N132D shows a clear shell feature from 
which most of the ion emission lines arise and a centrally confined Fe K emission charateristic 
of a high temperature component (Behar et al. 2001; Aschenbach 2002). According to the 
ASCA data analysis (Hughes et al. 1998), the high X-ray luminosity (~ 3 x 10 37 ergss _1 ) of 
the shell comes from the shock interaction with the dense material at the cavity wall, with 
an X-ray spectral normalization factor N s ~ 5 x 10 12 cm -5 , and hence a volume emission 
measure EM = And 2 N s ~ 1.5 x 10 60 cm~ 3 (where d « 50kpc) (Hughes et al. 1998). The 
ASCA spectral fit gives a temperature for the X-ray emitting gas kTx ~ 0.7 keV. Recently 
the XMM RGS spectrum indicates a low temperature component of ~ 0.6 keV (Aschenbach 
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2002), and the XMM EPIC-PN spectrum indicates ~ 0.9 keV (Behar et al. 2001). Compared 
with the central high temperature Fe K emission, the low temperature component could 
basically ascribed to the shell gas. If 0.7 keV is taken as a mean value of the postshock 
temperature, it corresponds to a shell velocity v s ~ 770 km s -1 , consistent with the above 
value derived from the optical observation. 

We consider a scenario in which the blast wave propagates rapidly in the wind cavity 
until it encounters the cavity wall, where it is slowed down drastically. Assuming the swept- 
up cavity material has been compressed sufficiently, the volume emission measure, under the 
thin-shell approximation, is given by 

EM « ^7rr s 3 [l - Aj(l - /?)] , (36) 

where the radius of the remnant r s is about 12 pc. For various values of f3 and Aj, we get 
from EM the estimates of the density at the cavity wall, n , and hence the swept-up gas 
mass, M s = 100M Q M 2 (eq. [1]), as plotted in Fig. 6 together with other physical parameters 
obtained from the following calculations. With the aid of eq.(8) for function i^ A) (Aj), we then 
get from eq.(7) an estimate of the explosion energy E . The age of the remnant, t = 10 3 t 3 yr, 
consists of two parts, tj (Sedov evolution before the blast wave reaches the cavity wall) and 
t—tj (duration after the impact). The latter is obtained from eq.(10) with eq.(15) for function 
^ (A) (Aj). 

It can be found from Fig. 6 that the kinematic age i ~ 3 x 10 3 yr and the preshock 
density n ~ 3 cm -3 can be reproduced by means of the wind cavity model, if parameters 

P ~ 0.1 and Aj ~ 0.94-0.95 

are adopted. With this set of parameters, we also have E 5 i ~ 3 and thus avoid an unreason- 
ably high inference. As a comparison, if f3 is adopted lower, E 5 i would be too high; and if Aj 
is larger, t would be too large compared with the kinematic age. If n is fixed at 3 cm -3 , the 
EM (eq.[36]) entails a density contrast (3 < 0.2. We note that, the previous overestimates 
of E corresponds to the case that i^ A) (Aj) is taken as 1 in eq.(7); in order for E to be of 
normal value, F^ A) (Aj) has to be > 1. This condition for F^ A) is satisfied around our best-fit 
parameters (3 ~ 0.1 and Aj ~ 0.94-0.95 (see Fig.l), which implies that the blast wave is still 
in the course of abrupt deceleration within the wind-blown shell. 

With the best-fit parameters, we also know that it takes the blast wave tj ~ 2.2 x 10 3 yr 
to travel in the cavity and that the velocity of the blast wave just before impacting the 
cavity wall is Vj ~ 1.9 x 10 3 kms _1 which is decreased to the present value ~ 800 km s" 1 
over a duration of ~ 700 yr. Also with these quantities, it is easy to confirm a sufficient 
compression of the swept-up cavity material, consistent with the above EM approximation 
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(eq.[36]). Moreover, the total swept-up mass is M s ~ 17OM , to which the gas originally 
inside the cavity contributes ~ 6OM . The latter mass seems too high for a wind cavity, 
however this is not a surprise, if a small part of gas are evaporated from the cavity wall or/and 
if interstellar clumps have been left within the cavity so that the cavity might contain a 
substantial amount of interstellar material. The existence of interstellar clumps is reasonable 
in view of the proximity of the remnant to the southern molecular cloud (Banas et al. 1997). 
In fact, a shocked interstellar cloudlet has been found in the remnant by the HST (Blair et 
al. 2000). In addition, the swept-up mass of the cavity gas (~ 6OM ) suggests that the SNR 
has been in the Sedov phase by the time it collides with the cavity wall. 

5. Conclusion 

Using the thin-shell approximation, we have developed a semi-analytic model for su- 
pernova remnants which evolve crossing a density jump in the environmental medium. The 
generic evolutionary relations are found for two cases, impact of the blast wave on a cavity 
wall and breakout of the blast wave from a dense cloud. In the impact case, both the adi- 
abatic and radiative phases are investigated; while in the breakout case, only the adiabatic 
phase is considered and the evolution relations are found to be an extension of those in the 
former case with different density contrast. In the impact case, it is also found that the 
remnant will evolve rapidly into the radiative phase and even the adiabatic phase could last 
so short that it seems to be abortive if the medium density at the cavity wall is sufficiently 
high. The developed model is applied to the cavity-born supernova remnant N132D whose 
evolution has not yet been quantitively estimated in a cavity scenario due to lack of proper 
model formulae and self-consistent physical parameters are obtained. 
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Fig. 1.— Plots of functions i^ A) (Aj) and F r (A) (Aj) for the impact model {(3 < 1). The solid 
lines denote the accurate values of eqs.(8) and (11). The four dotted lines, from upper to 
lower, stand for eq.(12) and eq.(16) in the limit case f3 = 0, and eq.(15) for (3 = 0.1,0.3, 
respectively, while the two short dashed lines stand for eq.(14) in the cases (3 = and (3 = 0.3. 
The long dashed line is the plot of eq.(18). 
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Fig. 2. — Plots of functions F^ (R) (A C ) and F r (R) (A c ) for rj — 1. The solid lines denote the 
accurate values of eqs.(27) and (30). The dotted lines stand for eqs.(31) and (35) in the case 
(3 = and eq.(34) in the cases (3 = 0.1,0.3, respectively. The dashed lines stand for eq.(32) 
for (3 = and eq.(33) for (3 = 0, 0.1, 0.3, respectively. 
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Fig. 3. — Plots of functions i^ (R) (A c ) and F r (R) (A c ) for rj = 1.2. The solid lines denote the 
accurate values of eqs.(27) and (30) and the dashed lines represent eq.(33). The solid and 
dashed curves of F r (R) (A c ) are plotted for (3 = 0,0.1, and 0.3, respectively, from upper to 
lower. 
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Fig. 4. — Plots of functions i^ (R) (A c ) and F r (R) (A c ) for r] — 1.5. The solid lines denote the 
accurate values of eqs.(27) and (30) and the dashed lines represent eq.(33). The solid and 
dashed curves of F r (R) (A c ) are plotted for (3 = 0,0.1, and 0.3, respectively, from upper to 
lower. 
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Fig. 5. — Plots of functions i^ R) (A c ) and F r (R) (A c ) for the breakout case {(3 > 1). The 
solid lines denote the accurate values of eqs.(8) and (11). The dotted lines represent the 
approximate series eq.(15). The long dashed line of F r (A) for f3 — 1 is plotted merely for 
comparison. 



-20- 




Fig. 6. — Plots of physical parameters derived for SNR N132D for (3 = 0.1 and Aj = 0.94, 
respectively. 



